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We study the electronic and transport properties of a graphene-based superlattice theoretically by using an 
effective Dirac equation. The superlattice consists of a periodic potential applied on a single-layer graphene 
deposited on a substrate that opens an energy gap of 2A in its electronic structure. We find that extra Dirac 
points appear in the electronic band structure under certain conditions, so it is possible to close the gap 
between the conduction and valence minibands. We show that the energy gap Eg can be tuned in the range 
0 < Eg < 2A by changing the periodic potential. We analyze the low energy electronic structure around 
the contact points and find that the effective Fermi velocity in very anisotropic and depends on the energy 
gap. We show that the extra Dirac points obtained here behave differently compared to previously studied 
systems. 


I. INTRODUCTION 

Graphene has attracted a great deal of attention since 
its first successful experimental fabricatioi>i in 2004 due 
to its intriguing physics and application potential^^— . 
Graphene is a one-atom thick layer of carbon atoms ar¬ 
ranged in a hexagonal structure and its low-energy elec¬ 
tronic structure can be described by using a Dirac-type 
Hamiltonian. The neutral, clean system has no gap and 
it is described by a massless Dirac equation. Due to the 
Klein tunneling^i^, charge carriers can not be confined by 
electrostatic potentials, what limits the uses of graphene 
in electronic devices. Opening a gap in the spectrum can 
help to confine the charges. 

An energy gap can be induced in graphene, for in¬ 
stance, by doping with boron^i^ or nitrogen^ atoms. An¬ 
other way to open an energy gap in the electronic struc¬ 
ture of graphene is using an appropriate substrate. It 
was verified that a hexagonal boron nitride (/i-BN) sub¬ 
strate induces an energy gap of 53 meV in graphenei^, 
which can be tuned by transverse electric fieldii. Epi¬ 
taxial graphene grown on SiG substrate has a gap of 
Ri 0.26 eV— . The other electronic property of graphene 
that depends on substrate is the Fermi velocity!^. The 
Klein tunneling can be suppressed also by electromag¬ 
netic fieldsi^^— and by a spatially modulated gapi^r— , 
leading to confined states. 

In the last years, the possibility of engineering the elec¬ 
tronic band structure of graphene by applying a peri¬ 
odic potential, i.e., a superlattice, has attracted consid¬ 
erable research interest to this subject. There are dif¬ 
ferent methods to generate the periodic potential struc¬ 
ture in graphene, such as electrostatic potentials^^Ti^S and 
magnetic barriers^^— . The combined effects of elec¬ 
trostatic and magnetic barriers have been studied as 
wel l^^i3^ . Despite the difficulty of fabricating graphene 
under nanoscale periodic potentials, it was already real¬ 
ized experimentally^i"— . It was found that the periodic 


potential leads to the appearance of extra Dirac points in 
the electronic structure of graphene^iiS^— and af¬ 
fect the transport properties, inducing an anisotropy in 
the carriers group velocitj^i^, leading to the collima- 
tion of electrons beam o^3i45,46 ^ electronic structure 
of a bilayer and trilayer graphene superlattice were also 
analyzed^. Periodic potential can not open an energy 
gap in graphene. 

In this paper, we investigate the electronic and trans¬ 
port properties of a graphene sheet deposited on a sub¬ 
strate that opens an energy gap in its electronic structure. 
On top of it we apply an external periodic potential. Our 
work is centered into analyzing the electronic structure 
in the vicinity of the new Dirac points that arise by the 
interplay of the gap induced by the substrate and the 
driven periodic potential. We show that the gap can be 
tuned by the external periodic potential and that the new 
Dirac points show characteristic differences with respect 
to those found previously in similar systems. We analyze 
the electronic and transport properties in the vicinity of 
the contact points be obtaining the dispersion relation 
and the effective Fermi velocity, which turns out to be 
very anisotropic around the contact points and is sensi¬ 
tive to the energy gap. 

The paper is organized as follows: In Sec. II we obtain 
the dispersion relation for the gapped graphene with a 
piecewise constant periodic potential. In Sec. Ill we in¬ 
vestigate the electronic and transport properties of the 
system. We analyze the electronic band structure for 
equal and unequal well and barrier widths and investi¬ 
gate the emergency of extra Dirac points. We also find 
the dispersion relation and the group velocity around the 
contact points. The paper is summarized and concluded 
in Sec. IV. 


II. THE DISPERSION RELATION 
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The electronic structure of a graphene sheet in the 
vicinity of a Dirac point K can be described by an ef¬ 
fective Dirac Hamiltonian. Applying an external one- 
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dimensional square-wave potential V (x) and considering 
an energy gap 2A in the electronic structure of graphene, 
the Dirac-like Hamiltonian reads 

ff = -ihVF{(7xdx + O-ydy) + V{x)l + AtJz , (1) 


Thus, we can define new effective mass and effective ky 
terms 

J A* = Acoscf) — hvpky sin (j) , , 

fc* = Asin^-I-cos(/> ' ' 


where Ui are the Pauli matrices, 1 is the 2x2 unitary ma¬ 
trix and uf is the Fermi velocity. The energy gap can be 
realized taking advantage of the influence of the substrate 
on the electronic properts of graphene. One can, for in¬ 
stance, deposit the graphene sheet in a SiC substratei^, 
as shown in Fig. O- Other important electronic prop¬ 
erty of graphene that is affected by the substrate is the 
Fermi velocity vp- For graphene in different substrates 
the Fermi velocity has been measured by different authors 
and their results summarized inH. We are considering a 
periodic potential V(x) with period a + b that is equal to 
V at 0 < X < a and zero at a < x < a -I- 6. 



FIG. 1. The graphene sheet deposited on a SiC substrate with 
an external periodic potential V (x) that has a period of a -|- 6. 
The substrate opens a gap of 2A in the electronic structure 
of graphene. 

The Dirac equation is given by 


Now, defining A/ihvpky) = tan(/) we make A* = 0. So, 
the Hamiltonian ([4]) is reduced to 

H'=-ihvpcfxdx + k^Uy+ V{x)i , (8) 

which is the two-dimensional massless Dirac equation for 
graphene with a periodic potential. The equation 


H'tj;'{x,y) = Eijj'{x,y) (9) 

was already solved by different method o^^i"*^^ and the dis¬ 
persion relation is given by 

cos{kxl) = cos{k*a) cos{k 2 b) 

(k*f+E(V-E) 

+ n'^y^k*k* -sm(fc*a)sin(fc*6) , (10) 

where fej) = {[E"^ — k^ = ([(H — E)"^ — 

{k*yf]/b?vlf/‘^, kx is the Bloch wave number and we 
have defined I = a + b. In order to transform back to the 
original ky and A terms one can use the inverse transfor¬ 
mation 


f A = A* cost/) -I- fc* sint/) 
yhvpky = — A* sint/)-b A:* cost/) 


( 11 ) 


As we define A/{hvpky) = tant/i, we have that {ky)‘^ = 
A^ + ti^Vpky. Replacing this in Eq. (nni we obtain that 
the dispersion relation for a 2D massive Dirac equation 
with a periodic potential is given by 


cos{kxl) = cos{kia) cos{k 2 b) 
kltv^vjr + E{V -E) + A^ 

h?Vpkik2 


sin(fcia) sin(A: 2 f)) (12) 


Htp{x,y) = Eip{x,y) , (2) 

where '0(x,j/) is a two-component spinor that represents 
the two graphene sublattices. Writing 

tp{x,y) = e*^“*'^/)(x) (3) 

and replacing in ([2]), one will have 

H = —ihvpOxd'x + hvpkyOy + H(x)l -I- Actz . (4) 

Applying the unitary transformation Ux{4>/‘^) = 
which commutes with ax and 1 but not with ay^z on can 
write 

H' = UxHUl = -ihvFCTxdx 

-1-6*^'^“’ [hvFkyUy + Acrz]e“*^ +V{x)i (5) 
Using the property Ujai = ii i ^ j we obtain that 

Ei'^^[hvFkyay + Aaz]e~"^'^^ = 

[A cos (j) — hvpky sin (/)] Gz + [hvFky cos (j) + A sin (/)]a^6) 


where fci = ([U^ — A'^]/h^Vp — ky)^^'^ and k 2 = ([(U — 
EY — A’^]/b?v‘]^ — ky)^/"^. Note that at U = A = 0 we 
recover the linear dispersion relation of a graphene sheet. 

The left hand side of Eq. (TT^ is limited to the interval 
(-1,1). Therefore, in the right hand side one has allowed 
and forbidden values for the energy, which implies in the 
appearance of energy bands with gaps. 

III. THE ELECTRONIC STRUCTURE 

Having obtained the dispersion relation, in this sec¬ 
tion we will analyze the electronic structure. In what 
follows, we shall consider a constant period of the super¬ 
lattice equal to 60 nm, i.e., a -I- 6 = 60 nm. As Eq. m 
is invariant under simultaneous replacements E —>■ —E 
and U —>■ — U, only non-negative values of V will be 
considered. We shall concentrate our discussion on the 
valence and conductance minibands only, assuming the 
Fermi level to be in between at any value of V. 
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A. Band structure with equal well and barrier widths 

Here we will study the electronic band structure at 
a = 6 = 30 nm in the dispersion relation m, which 
means that the well and barrier have the same width. In 
Fig. [2] (a) are plotted the electron and hole energies as a 
function of with ky = 0 and A = 0.13 eV for V = 270 
meV (red), V = 304.52 meV (black) and V = 340 meV 
(blue). It can be seen that when the potential increases, 
the electron and hole minibands shift up. However, the 
shift of the electron miniband is not equal to the shift of 
the hole miniband, which implies different electron-hole 
minigaps for different values of H, as shown in Fig. [5] 
(a). One can see that it is possible to close the minigap, 
as happens when V = 304.52 meV (black), showing that 
is possible to have a gapped or gapless graphene only 
changing V. It is a consequence of having a position 
dependent potential. If the potential is constant, the 
electron and hole minibands are shifted equally and the 
minigap remain the same, regardless of the value of the 
potential. 




y(eV) 


FIG. 2. The dispersion relation (US with a = b — 30 nm. (a) 
The electron and hole energies as a function of kx , with ky = 0 
and A = 0.13 eV at F = 270 meV (red), V = 304.52 meV 
(black) and V = 340 meV (blue), revealing the possibility of 
closing the energy gap with the potential, (b) The electron 
and hole minibands as a function of V with kx = ky = 0, 
which shows the oscillation of the energy gap. 

This is more clear when we look to Fig. [5] (6), where 
the electron and hole energies are plotted as a function of 
V with kx = ky = 0. It should be noted that the contact 


point is obtained in kx = 0, as can be seen in Fig. 
(a). Therefore, Fig. [5] (6) is showing how the electron- 
hole minigap changes with the potential for ky = 0. It 
can be seen that the minigap oscillates when V changes, 
and may be zero. For the values of the parameters chosen 
here, the first value of V that closes the gap is F = 304.52 
meV. The highest value for the energy gap is obtained 
at F = 0, which is equal to 2A. Thus, it is not possible 
to increase the initial gap in graphene with a periodic 
potential. It means that a periodic potential can tune 
the Dirac gap Eg only in the range 0 < Eg < 2A. So, if 
A = 0, the potential is not able to open a gap. 



FIG. 3. The electron and hole energies as a function of ky 
with a = 6 = 30 nm, kx ~ 0 and A = 0.13 eV at F = 270 
meV (red), F = 304.52 meV (black) and F = 340 meV (blue). 
Extra Dirac points appear at fcj, 7 ^ 0 when F exceeds a critical 
value Vc- In this case Vc = 304.52 meV, and the energy gap 
remains closed. 


The electron and hole energies as a function of ky with 
kx = 0 are plotted in Fig. [3] for the same values of V as 
in Fig. [2] (a). When V = 270 meV (red) and V = 304.52 
meV (black) the minibands have the same behavior that 
in Fig. (a), but are narrower. For V = 340 meV 
(blue) the minigap at ky = 0 opens, however there are 
extra Dirac points appearing at different values of ky. 
These extra Dirac points appear when V exceeds a crit¬ 
ical value, that is Vc = 304.52 for the values of the pa¬ 
rameters chosen here, and do not disappear. So, from 
V = Fc, the gapped graphene becomes gapless. 

In order to find an expression for Vc in terms of the 
system parameters, let us first localize the contact points 
in k space. Taking into account the implicit function 
theorem, one can conclude that at the contact points, 
where there is an intersection of the bands, the gradient 
(Jacobian) of the dispersion relation should be zero. Note 
that ki = k^ when E = Eq = V/2 and that the contact 
points are all at kx = 0. So, the Eq. (HU with a = b, 
kx = 0 and E = Eq is given by 


I = cos^(A:ia) -I- 


k^h^vl + Eo{V - Eq) + 


fi^Vpkl 


sin^(fcia) 


(13) 

which is satisfied when ki = mr/a, where n is an integer 
different of zero, because n = 0 implies fci = 0, which 
makes the denominator b?Vpkf in Eq. (US vanishes. 
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This condition leads to 



which gives the values of ky where the contact points 
are located. The exact location of the contact points 
are {E,kx,ky) = {Eo,0,ky^). Should be remembered 
that the Dirac points appear only after a critical value of 
V = Vc. The zeros of the equation above give the contact 
points at ky = 0. So, one can write 


14 = 2i 


I ^ mrhv ^' 


+ A2 


(15) 


which gives the values of V where there is a contact point 
in Fig. [5] (b). The critical potential is given by 14 = 14- 
The number of contact points can be found from Eq. 
(HI- When {Eq — A'^y^'^a/hvp'^ is not an integer, the 
number of contact points is given by 


Nd 


a^jEl - A2 
Trhvf 


(16) 


where [• • •] denotes an integer part. When {Eq — 
= {mr/a)^, the number of Dirac points is 
Njj = 2n — 1. A different way to obtain the number 
of Dirac points is: when Vn < V < 14+i, Ejj = 2n, 
whereas when V = Vn, Njj = 2n — 1. 


B. Band structure with unequal well and barrier widths 




FIG. 4. The dispersion relation (HU with a ^ b. The dashed 
lines are the minibands with a = 20 nm and & = 40 nm, 
whereas the continuum lines are the minibands with a = 40 
nm and 6 = 20 nm. (a) The energy in terms of kx with ky = Q 
and A = 0.13 eV at F = 290 meV (red), V = 319.07 meV 
(black) and V = 340 meV (blue). (6) The electron and hole 
minibands as a function of V with kx = ky = 0. 


Now let us consider the case with a ^ b. The dashed 
lines in Fig. 3] are the minibands with a = 20 nm and 
6 = 40 nm, whereas the continuum lines represent the 
minibands with a = 40 nm and 6 = 20 nm. In Fig. |4l(a) 
are plotted the electron and hole minibands as a function 
of kx with ky = Q atV = 290meV (red), V = 319.07meV 
(black) and V = 340meV (blue). As in the case with a = 
b, for different values of V there are different electron- 
hole minigaps, which may be zero. The oscillation of the 
minigap is shown in Fig. |4] (5), where the minibands as 
a function of V with kx = ky = 0 are plotted. 

It can be seen that the values of V that close the gap 
when a = Cl and b = C 2 are the same for a = C 2 and 
b = Cl- However, due to the fact that the potential shifts 
the minibands, when the graphene region with V{x) = V 
is wider than the region with V{x) = 0, there is a larger 
shift of the minibands. It explains the difference in energy 
between the dashed and continuum lines in Fig. |4l 

In Fig. [S] we plotted the minibands as a function of 
ky with kx = 0- Once more, the dashed lines are the 
minibands with a = 20 nm and 6 = 40 nm, whereas the 
continuum lines are the minibands with a = 40 nm and 
& = 20 nm. In Fig. [S] (a) we recovered the first time that 
the minigap closes at F = 319.07 meV. In Fig. [S] (6) we 
have V = 400 meV. One can see that, when a ^ b, the 


extra contact points that appears at ky 0 are not in 
the Fermi level. When a > b {a < b) the contact points 
are shifted up (down) the Fermi level. 

In order to localize the contact points, again, we take 
advantage of the implicit function theorem. The gradient 
of the dispersion relation will be zero only if sinfcio = 
sinfc 2 & = 0 and cosfcio = cosfc 2 & = if- So, one can 
write 


kia 


- A2 

V 



(17) 


and 


k2b = 


({V -Ef 



b = mn , 


(18) 


where m is an integer. Subtracting (1181) from dm , one 
gets 


F ■n’^b^v^p / w? rn?\ 

^ Y 2F ~1^) 


(19) 


Replacing the equation above in Eq. (II3 one obtains 


= ±A 


fl^Vp 


( 20 ) 
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FIG. 5. The electron and hole minibands as a function of ky 
with A = 0.13 eV and fca; = 0 at (a) V = 319.07 meV and 
(6) V = 400 meV. The dashed lines are the minibands with 
a = 20 nm and 6 = 40 nm, whereas the continuum lines are 
the minibands with a = 40 nm and 6 = 20 nm. The extra 
Dirac points at 7 ^ 0 are shifted up (down) the Fermi level 
when a > b {a < b). 

From the zeros of equation above one obtains, 

which is the values of V where there is a contact point 
in Fig. 0] b. Again, the critical potential, where the 
graphene superlattice becomes gapless, is 14 = Vi. One 
can see that increasing the difference between a and b the 
value of 14 increases, as well. Therefore, for a particular 
V, there is always a value of a and b which the graphene 
is gapped. When a = 6 in Eqs. dUll-lEIl), the results 
obtained in the last section are recovered. 

The exact location of the contact points when V = Vm 
OT Vm <V < Vm+i is given by 

{E, ky) = {Em, ky^), {Em-1, {^1, ^Vl) , (22) 

where {Em,ky^) is the contact point nearest to ky = 0, 
whereas {Ei, ky^) is the contact point farthest to ky = 0. 
Remember that in our system all contact points are at 
kx = 0- An example is given in Fig. [6l where the electron 
and hole energies are plotted as a function of ky with 
a = 20 nm, 6 = 40 nm and V = V 3 = 600.1 meV. 
The location of the contact point nearest to ky = 0 is 



FIG. 6. The five Dirac points in the k space with a = 20 nm, 
& = 40 nm and V = V 3 = 600.1 meV. The Dirac points at 
ky ^ Q are located below the Fermi energy. 


(i?3,0). The next nearest contact point is {E 2 ,ky 2 ) and 
the farthest contact point is located at {Ei,ky^). 

The number of contact points is the same obtained in 
the last section with a = b. 


C. Dispersion relation near the contact points 


Now, let us analyze the electronic structure near the 
contact points. The electronic structure in the vicinity 
of the contact points has been studied in the case of a 
gapless graphene superlattice^^ and for a graphene su¬ 
perlattice with spatially modulated gap^. In both cases, 
the discussion was restricted to the particular case a = b. 
For the sake of comparison, we will consider the special 
case with a = b and then consider the general case with 
a ^ b. For this purpose, one has to expand Eq. (fT^ 
in the vicinity of the contact points obtained in the last 
sections. 

Let us first consider the special case a = b. In order 
to obtain the behavior of all contact points, one has to 
consider the contact points at E = Eo{Vn), kx — 0 and 
ky = 0 and the contact points at E = Eo{Vn), kx = 0 
and ky = ky^. Thus, expanding Eq. (TT^ into the Taylor 
series up to second order of E — Eo{Vn), kx and ky, one 
gets 


ei = ± 


kp^d^u^ - 4 ^ 4 ) 
kyE{2u^ — 8d‘^) + 


— 256(i® -b 4it® — 48u^(i^) 
kyE{2u^ — 8d‘^) + — 167r^d^u‘^ 


1/2 

(23) 


where we have defined ei = {E — EQ)l/hvF, u = 
Vnl/A ttHvf and d = Al/AirhvF- The positive and neg¬ 
ative signs represent the electron and hole minibands, 
respectively. Writing ky = ky — ky^ and expanding Eq. 
m up to the lowest order of E — Eq, kx and ky, one 
obtains the dispersion law in the vicinity of the contact 
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points located at ky ^ 0 , which is given by 


e2 = ±WA:2F^ + fcJMl-2^+4 


, (P-vi?- 


-4 


d2 


(24) 

It should be noted that, when a = b, the electron and 
hole minibands are symmetric related to = 0 , which 
does not happen when there is a periodic modulation of 
the energy gap in graphene^. 




FIG. 7. The electron and hole energies in terms of ky with 
fca; = 0 in order to compare the exact dispersion relation (blue 
line) with the spectrum obtained around the contact points 
(red dashed line), (a) £i at V = Vi. (b) £2 at F = 340 meV 
and n = 1. (c) £3 at V = 600.1 meV and m = 3. (d) £4 at 
V = 600.1 meV and m = 2. 


Considering now the most general case with a 7 ^ 6 , one 
can expand Eq. (HU with V = Vm up to the lowest order 
of E — Em, kx and ky, and obtain the dispersion relation 
in the vicinity of the contact points at kx = ky = 0 and 
E = Em, that is given by 

£3 = aikl ± ^JPik^ + likl . ( 26 ) 

The dispersion near the contact points at ky ^ 0 can be 
written as 

£4 = a 2 ky ± \J /32kl + 72^2 , ( 27 ) 

where ky = ky — ky^. The coefficients ai, /3i and 7 ^, with 
j = 1, 2, depend on a, b, V and A. They are too large to 
be write down here. When a = b, the coefficients ai and 
/ 3 i vanish and £3 becomes iLkx(V,^ — 4A^)/I^. For this 
reason, it was necessary to expand Eq. m up to second 
order of E — Eo{Vn), kx and ky to get £ 4 . 

One can note that, in contrast to the case with equal 
well and barrier widths, when a ^ b the electron and 
hole minibands are not symmetric related to £i = 0 due 
to the coefficient ai. However, as in the case with a = b, 
the minibands along the ky direction are parabolic in the 
contact points at ky = 0 and conical at ky 7 ^ 0. In Fig. [7] 
(a) and (b) we compare £3 and £4, respectively, with the 
exact dispersion relation We consider V = V3 = 

600.1 meV and plotted £4 in the vicinity of the contact 
points located at {E 2 ,ky.^). 

Should be mentioned that the dispersion relation is 
linear along kx around all contact points. So, the energy 
surface £{kx,ky) is conical in the vicinity of the contact 
points at ky ^ 0 and has a lenslike shape around the 
contact points at ky = 0. 


D. Group velocity around the contact points 


In Fig. [7] (a) we plotted £1 (dashed red line) with 
V = Vi and kx = 0 as a function of ky and compare with 
the exact dispersion relation ([T2l) (blue line). One can 
see that the expansion (1231) is good in the vicinity of the 
contact point. Expanding £1 in powers of ky with kx = 0, 
one gets 

> ( 25 ) 

which gives a parabolic electron and hole minibands, as 
can be seen in Fig. [7] (a), in contrast to the conical 
dispersion around the original Dirac point in a gapless 
graphene. In the limit when ^ d, the dispersion 
along ky becomes flat. Note that in the case of a gapless 
graphene superlattice^^, the dispersion along ky is given 
by £1 ifcy. In Fig. [7] (5) we compare £2 (dashed red 
line) with Eq. (ITT!) (blue line) at V = 340 meV. Again, 
there is an agreement between the exact dispersion re¬ 
lation and the expansion (1241) near the contact point. 
However, the dispersion is linear along ky, which does 
not happen in the contact points located at ky = 0. 


Let us now use the spectrum for small energies ob¬ 
tained above to find the effective Fermi velocity around 
the contact points. The components of the velocity in 
the vicinity of the contact points are given by Vx^vp = 
dsi/dkx and Vyjvp = dsi/dky, where i = 1,2, 3,4 de¬ 
note the four kinds of contact points. The expressions 
for the components of the velocity can be seen in the 
Appendix. 

The anisotropy of the electron and hole velocities in the 
(kx, ky) plane can be seen clearly if one introduces a polar 
angle ip with the relations kx = Q cos ip and ky = Q sin ip, 

where Q = ^k^ -(- ky. In Fig. | 8 ]the absolute value of the 

velocity Vi = as a function of ip was plotted 

for two different values of the energy gap: A = 0.13 
eV (continuum line) and A = 26.5 meV (dashed line), 
which correspond to the graphene on a SiC and /i-BN 
substrate, respectively. It should be mentioned that the 
Fermi velocity Ui? in graphene on these two substrate is 
1.15 • 10® m/s and 1.49 • 10® m/s, respectively^^. One 
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can see that the velocity is sensitive to the energy gap A. 
The velocity has smaller values for larger values of A. 



FIG. 8. The electron (blue) and hole (red) group velocities 
in the vicinity of the four kinds of contact points as a func¬ 
tion of (p at A = 0.13 eV (continuum line) and A = 26.5 
meV (dashed line), (a) and (b) vi/vp at Vcont = 304.52 meV 
and Vdash = 212.03 meV. (e) vij^F at Vcont = 357 meV and 
Vdash = 313 meV. In these cases, the electron and hole veloc¬ 
ities are equal, (c) and (d) vs/vf with a = 20 and 6 = 40 nm 
at Vcont = 319.07 meV and Vdash = 237.65 meV. (/) VifvF 
with a = 20 nm and 6 = 40 nm at Vcont = 383.5 meV and 
Vdash = 351.48 meV. Vcont and Vdash are the value of V for 
the continuum and dashed lines, respectively. 

In Fig. [5] (a) and (e) we plotted the velocity for a = 5 
at V = Vi (Fi = 304.52 meV when A = 0.13 eV and 
Vi = 212.03 meV when A = 26.5 meV) and at the inter¬ 
mediate value between Vi and V 2 {V = 357 meV when 
A = 0.13 eV and V = 313 meV when A = 26.5 meV), 
respectively. When a = b the electron and hole veloci¬ 
ties are equal, in consequence of the symmetry between 
the electron and hole minibands. This does not happen 
in the case with a modulated energy gap^, where the 
electron and hole velocities are not the same. For small 
values of A, is close to vp for almost all values of the 
angle </?, having a narrow dip in the vicinity of ip = 'n{2^ 
as can be seen in Fig. [5] (a). This is due to the fact that 


Vx-i is much greater than for all angle tp except in 
the vicinity of 7r/2, where both Vx^ and Vy-^^ are small, in 
consequence of the lenslike shape of the energy surface. 
When A increase the velocity decrease and the value of 
vi become much smaller than vp. The dip remains at 
tp = 7r/2, but it become a little wider. A similar behav¬ 
ior was obtained in^. In Fig. [5] (5) there is a zoom of the 
dip region. Different of ei, the energy surface generated 
by £2 is conical, but it is not an isotropic cone, generating 
an anisotropy in the velocity. In this case, with a small 
energy gap, the electron and hole velocities have only a 
little variation around vpl2, as can be seen in Fig. [5](e). 
Increasing A, there is a stronger anisotropy, in contrast 
withes. 

The electron and hole velocities for a = 20 nm and 
5 = 40 nm at F = Fi (Fi = 319.07 meV when A = 0.13 
eV and Fi = 237.65 meV when A = 26.5 meV) and at 
the intermediate value between Fi and F 2 (F = 383.5 
meV when A = 0.13 eV and F = 351.48 meV when 
A = 26.5 meV) are plotted respectively in Fig. [H](c) and 
(/). When a ^ b the electrons and hole minibands are 
asymmetric, so the electron and hole velocities are not 
equal. In Fig. [5] (c) we consider the contact point at 
ky = 0. The behavior of the velocity in this case is very 
similar with the case with a = b. The main difference 
is that the dip has a width slightly different. One can 
note that the electron and hole velocities are almost the 
same, differing slightly in the vicinity oi ip = p 12. The 
Fig. [5] (d) is an extension of the dip region. The energy 
surface generated by £4 is conical, but is a tilted and not 
isotropic cone. In Fig. [5] (/) we plotted V4. It can be 
seen that the electron and hole velocities are equal at 
ip = 0,p and differ widely for other values of ip. As in 
the case with a = b, when A increases the anisotropy of 
the velocity becomes greater. 


IV. CONCLUSIONS 

We have analyzed the electronic structure of a gapped 
graphene superlattice with a piecewise constant periodic 
potential using the continuum model based on an effec¬ 
tive Dirac equation. We consider that the energy gap 
is generated by an appropriate substrate, which changes 
the Fermi velocity, as well. 

It was shown that the energy gap oscillates when the 
potential F changes continuously at kx = ky = 0 and is 
zero at discrete values Fi- When F > Fi, extra Dirac 
points appear at fcj, 7 ^ 0 and never disappear. Thus, be¬ 
ginning with a critical potential Fc = Fi, the graphene 
system becomes gapless. In the special case of equal 
well and barrier widths, these extra Dirac points are lo¬ 
cated in the Fermi level and the electron and hole mini¬ 
bands are symmetric, whereas with an unequal well and 
barrier widths the extra Dirac points are no longer at 
the Fermi level and the electron and hole minibands are 
asymmetric. We found that if the initial energy gap Eg in 
graphene is equal to 2A, it is possible to tune the energy 

























gap with a periodic potential in the range 0 < Eg < 2A. 
We found the locations of all contact points and it was 
shown that the greater the difference between the well 
and barrier width, the greater the critical potential 14- 
Finally, we obtained the dispersion relation in the vicin¬ 
ity of all contact points and used it to find the effective 
group velocity of the carriers. The velocity has a strong 
anisotropy around the contact points and is sensitive to 
the energy gap. Extra Dirac points of different kinds have 
been already studied in graphene superlattices. Analyz¬ 
ing the electronic structure near the contact points, we 
showed that the extra Dirac points obtained here have a 
different behavior compared to previously studied. The 
results obtained here can be used in the fabrication of 
graphene-based devices. 



8 

(A.7) 

(A.8) 

(A.9) 

(A.IO) 


ACKNOWLEDGMENTS 

I thank M. A. H. Vozmediano for helping me to re¬ 
vise and correct mistakes in a previous version of the 
manuscript. This work was partially supported by CNPq 
and CNPq-MICINN binational. 

Appendix: The components of the group velocity around 
the contact points 


The components of the velocity in the vicinity of 
the contact points are given by v^Jvp = detjdkx and 
Vyilvp = dsi/dky, with i = 1,2,3,4. Thus, 



C = l‘^{2u^ (A.ll) 

and 

D = — A(f) . (A-12) 
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